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^ ' Abstract 

^ We point out that the existence of post-inflationary phases stiffer than radiation leads to the 

^ , production of stochastic gravitational waves (GW) backgrounds whose logarithmic energy spec- 

Q^ , tra (in critical units) are typically "blue" at high frequencies. The maximal spectral slope (for 

CO ' present frequencies larger than 10~^^ Hz) is of order one and it is related to the maximal sound 

Q velocity of the stiff plasma governing the evolution of the geometry. The duration of the stiff 
00 

phase is crucially determined by the backreaction of the GW leaving the horizon during the 
^ de Sitter phase and re-entering during the stiff phase. Therefore, the maximal (inflationary) 
I curvature scale has to be flne-tuned to a value smaller than the limits set by the large scale 

_fi y ^ 

O , measurements (ifrf5<10 Mp) in order to a have a sufficiently long stiff phase reaching an 
■ energy scale of the order of the 1 TeV and even lower if we want the stiff phase to touch the 
^ ' hadronic era (corresponding to Thad ~ 140 MeV). By looking more positively at our exercise we 
- see that, if an inflationary phase is followed by a stiff phase, there exist an appealing possibility 
of "graviton reheating" whose effective temperature can be generally quite low. 



^Electronic address: m.giovannini@damtp.cam.ac.uk 



1 Introduction 



Strong causality arguments forbid the existence of a never ending radiation dominated epoch. 
If this would be the case regions emitting a highly homogenous and isotropic Cosmic Microwave 
Background Radiation (CMBR) at the decoupling epoch, would not have been in causal contact 
in the far past. This problem of horizons (together with other kinematical problems of the 
standard cosmological model) led to assume an inflationary phase of accelerated expansion 
(a > 0, a > 0) where the (effective) equation of state describing the background sources during 
inflation was drastically deviating from the one of radiation. 

In the context of ordinary inflationary models the transition from the inflationary regime 
to the radiation era is usually associated with a reheating phase where the energy density of 
the inflaton field is released, "producing" a radiation dominated phase. The dynamics of the 
inflaton right after the inflationary epoch has been recently discussed in detail 0. 

If this if the dynamical picture of the evolution of our Universe in its early stages, then, 
one of the most interesting (and difficult to test) implications is the production of a stochastic 
background of gravitational waves (GW). It is indeed known since many years that the various 
transitions of the curvature scale lead, necessarily, to the amplifications of the quantum me- 
chanical (vacuum) fluctuations of the tensor modes of the geometry |Q and to the consequent 
production of highly coorelated graviton squeezed states |^ . According to this mechanism, 
the energy density of gravitational origin can be estimated to be (today) of the order of 10~^^ 
[^, ^, 3 (in critical units) for (present) frequencies larger than 10~^^ Hz. The (logarithmic) 
energy spectrum turns out to be flat in the same interval of frequencies. The reason of this 
quite minute amplitude comes essentially from the measurement of the CMBR anisotropics. In 
fact, the tensor contribution to the CMBR anisotropics imposes a quite important bound on 
the maximal curvature scale at which the inflationary expansion occurred. By assuming, for 
example, that the inflationary phase was of de Sitter type, with typical curvature scale Has, 
we have to require [|1^, pJ]] , 

Mp^ ' ^ ' 

in order to be compatible with the detected level of anisotropics in the microwave sky. Con- 
straints on the dynamical assumptions of various inflationary models can be derived on the 
same basis |I2| . 

In this paper we want to explore a slightly different picture of the post-inflationary phase. 
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Our suggestion is, in short, the following. Suppose that the inflationary phase is not immedi- 
ately followed by a radiation dominated phase but by an intermediate phase whose equation of 
state is stiffer than radiation (i.e. p = jp, with 7 > 1 and 7 = [cs is the sound velocity of the 
plasma]). Then, two interesting implications can arise. On one hand the transition between 
the inflationary regime and the stiff regime leads to graviton spectra which slightly increase 
with frequency (with "blue" slopes) in the ultra-violet branch of the spectrum, on the other 
hand, the backreaction effects associated with the GW leaving the horizon during the de Sitter 
phase and re-entering during the stiff phase can heat up the Universe very efficiently, leading 
to a qualitatively new kind of reheating which one can call "graviton reheating" . 

Even if the implications of this suggestion can be appealing it is certainly important to 
justify better some possible motivations of such a weird exercise. 

First of all we can say that such a suggestion is not forbidden by any present data. An 
indirect evidence of the fact that the Universe might have been dominated by radiation around 
temperatures of the order of 0.1 MeV comes from the success of the simplest (homogeneous 



and isotropic) big-bang nucleosynthesis (BBN) scenario |[T3|. In the absence of any external 
magnetic field [111] and of matter-antimatter domains [|1^], the light elements (^ He, ^ He, Li, D) 



abundances are reproduced by the BBN model provided the ratio of the baryonic charge density 
over the photon density is fine tuned around 10^^'^. However, prior to nucleosynthesis, there 
are no direct tests of the thermodynamical state of the Universe and, therefore, the effective 
equation of state of the perfect fluid sources driving the evolution of the background geometry 
can be arbitrarily different from the one of a radiation dominated plasma. Moreover, there are 
no compelling reasons why long-range (Abelian) gauge fields should not have been present in 
the early Universe []T3[. Owing to our ignorance of the thermodynamical state of the Universe 
prior to the nucleosynthesis epoch, it is possible to postulate, in the framework of a particular 
model, the existence of post-inflationary (decelerated) phases different from radiation. 

The possibility of having post-inflationary phases whose (effective) equation of state was 
stiffer than radiation can be also motivated in terms of different models. It was recently argued 
that a stiff phase could be originated by the relaxation of the moduli towards the minimum 
of their non-perturbative potential [0. A similar idea was investigated in the framework of 
the dilaton relaxation |]1^ where the resonant amplification of gauge fields modes was also 
discussed. In fact one can argue that it does certainly exist a regime where the kinetic energy 
of the single modulus is dominant against its potential energy. Therefore, the geometry would 
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evolve following an effective equation of state whose sound velocity is approximately equal to 
the speed of light. In this regime the energy density of the modulus would scale as pmod ~ 
Of course intermediate situations can be also imagined so that we can, say, more generally that 
some phases with 1/3 < 7 < 1 can occur as a result of the moduli relaxation. For example, in 
the context of dilaton relaxation ||T^, if the kinetic energy dominates we get exactly a stiff fluid 
model with 7 = 1. In fact, on theoretical ground the dilaton field (0) has a potential going to 
zero (in the supersymmetric limit) as a double exponential (i.e. V ~ exp [— c^exp (—0)], with 
positive and model dependent constant). On more physical ground the V{(f)) is believed to 
have one (or more) minima for some ~ 0min- When reaches curvature scale if ~ m (where 
m is the dilaton mass) an oscillating phase begins. Depending upon the initial conditions in 
Hi (the maximal curvature) the oscillating phase can be preceded by a stiff phase where the 
dilaton kinetic energy dominates (0^ ^ V^(0)), and, therefore the evolution of the dilaton, in 
curvature, will be = 0o + 0i log [H/ Hi] implying 0^ ~ . 

It was recently speculated |TB[ (without relation to moduli relaxation) that, provided the 
stiff phase is long enough, interesting effects can also be expected at energy scales of the order 
of 100 GeV (corresponding to a curvature scale of the order H^.^ ~ 10^^^ Mp). Historically, 
the first one to imagine the appealing possibility of having long stiff epochs was Zeldovich |1^ . 
At that time inflationary models were not yet formulated and it seemed quite crucial to model 
correctly the quark-hadron phase transition occurring when the temperature of the Universe 
was of the order of the pion mass (i.e. T^ad ~ 140 MeV). The idea was that, at this stage, the 
equation of state of the perfect fluid sources could be stiffer than the one of radiation, namely 
with 7 > 1/3. On physical ground, one can easily understand that the sound velocity in a 
perfect fluid is likely to be smaller than the speed of light (7 < 1 in our units). In particular 
the stiff model of Zeldovich assumes exactly that, prior to the hadronic stage of evolution, 
7 = 1. The fact that the speed of light equals the speed of sound also implies that the energy 
density of the stiff sources decreases faster than radiation. 

This last feature of the stiff picture of hadronic interactions is indeed fatal for the logical 
consistency of the whole proposal. Actually, the existence of a hadronic phase with effective 
equation of state stiffer than radiation also implies the production of a stochastic background 
of GW sharply peaked towards the Planck frequency [^. Thus, the produced gravitational 
waves will backreact on the geometry (effectively driven by the stiff fluid). Now, since the 
energy density of the high frequency GW re-entering the horizon during the stiff epoch scales 



3 



like radiation ^ it was correctly concluded that, if a stiff phase ever existed prior to 
the usual radiation dominated phase, the back reaction effects associated with the production 
of high frequency gravitons turned very quickly the evolution of the Universe into a radiation 
dominated phase. It was also showed [^, that there are no chances of having a stiff 
phase from the Planck curvature scale {Hp ~ Mp) down to the hadronic curvature scale 
{H^a, ~ 10-40Mp). 

In our context there is a crucial difference with respect to Zeldovich suggestion and it is 
essentially given by Eq. ( [L . 1| ) . Since the curvature scale is quite minute in Planck units GW 
backreaction is not switched on immediately. None the less, the aim of this exercise is to point 
out that the occurrence (and the duration) of a post-inflationary phase stiffer than radiation 
is not a free parameter which we can adjust in the framework of a particular model to get 
the desired effects. On the contrary the duration of a stiff phase is significantly constrained 
by the back reaction of hard gravitons. In this sense our considerations owe very much to the 
pioneering works in the subject [^, |22[ but are applied to a different dynamical picture, where 
an inflationary phase is followed by a stiff phase. Our logic, in short, is the following. Let us 
assume, for example that in some specific model a de Sitter (inflationary) phase is followed 
by a stiff phase (for example with p = p) at some cosmic time ti. Then, the hard gravitons 
excited during the de Sitter phase will have typical (physical) momentum Tiujdsiti) ~ fi-Hds{ti)- 
Since the the Hubble distance deviates (during the stiff phase) from the value during the 
de Sitter epoch as a consequence of the change in the expansion rate the hard gravitons with 
^ ^ds will "re-enter" at different times during the stiff phase. Since their effective equation 
of state is the one of radiation they will modify the dynamics of the stiff phase by ultimately 
destroying it and by turning it into a radiation dominated phase. By looking positively at this 
effect we could say that the graviton back reaction represents a reasonable candidate in order to 
implement a reasonable reheating mechanism in these classes of models. At the same time the 
constraints imposed by the (over)-production of gravitons might be viewed as a weakness of the 
scenario. A related result of our analysis will be the calculation of the GW spectra produced 
during these types of stiff post-inflationary phases. 

The plan of our paper is then the following. In Sec. II we will introduce the basic equations 
describing the background evolution during a stiff, post-inflationary, phase. Sec. Ill is devoted 
to the calculation of the GW spectra in these models. In Sec. IV we will discuss the back 
reaction effects associated with hard (non-thermal) gravitons and we will discuss the possible 
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modification of tlie background evolution. Sec. IV contains our concluding remarks. 



2 Basic Equations 



In this Section we consider the simplest homogeneous and isotropic models of FRW type with 
line element 



ds = g^iydx^dx" = dr - a (t) 



+ r^(rfr + sin' 



(2.1) 



. 1 — nr^ 

(Greek indices run from to 3, whereas Latin indices run from 1 to 3). The derivative with 
respect to the cosmic time t will be denoted by an over dot whereas the derivative with respect 
to the conformal time t] will be denoted by a prime (as usual a{ri)dri = dt). The sign of the 
spatial curvature k = +1, 0, —1 corresponds to a closed, flat or open space. We will assume 
that the evolution of the geometry follows general relativity. Therefore the coupled evolution 
of the (perfect) fluid sources and of the geometry will be conveniently described in terms of the 
well known FRW equations 



Ml 



p + 3H{p + p) 



K 

-^(P + 3p), 



0, H 



d\oga 
dt '' 



(2.2) 



-1/2 



where H is the Hubble parameter and we also write Mp = |^87rG'/3j . Focusing our attention 
on the conformally fiat case (k = 0) we will consider models of background evolution where an 
inflationary phase (a > 0, a > 0) is followed at some time ti by a decelerated phase (a < 0, 
a > 0). In particular we will assume that the sources for t > ti will be well approximated by a 
barotropic equation of state 

P = 7P, (2.3) 

with 1/3 < 7 < 1. We consider unrealistic the case where 7 > 1 since this would mean that the 
sound velocity of the fluid is greater than the speed of light. By integrating the FRW equations 
( |2.2|) after ti with the closure given by Eq. (|2.3|) we get that 



a{t) = ai 



t 

ati 



P{t) = Pi 



ai 



(2.4) 
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with 

Pi = HlMl, ai = a{ti), (2.5) 

3(7 + 1) 

where Hi = H[ti) is simply the value of the Hubble parameter at the end of the inflationary 
phase. 

We want to stress that our approach in the present section will be an effective one: since the 
evolution equations of the tensor modes of the geometry is essentially determined only by the 
behavior of the scalar curvature we feel free of using the simplest fluid model for the description 
of the background sources. This is of course not a limitation. The effective fluid sources can 
be thought as modelled by the energy momentum tensor of one (or more) scalar fields. 

Concerning the inflationary phase we will not make any type of weird assumption about 
the background evolution. Indeed we will show that our considerations will be only mildly 
sensitive to the specific inflationary dynamics and will only (but crucially) depend upon the 
maximal curvature scale reached during inflation. From a purely kinematical point of view 
we will explore expanding inflationary epochs (i.e. d > 0, a > 0) with constant or decreasing 
curvature (i.e. H <0). The de Sitter (or quasi de Sitter) case seem to emerge quite naturally 
in the framework of the slow-rolling approximation 0. It is important to point out that the de 
Sitter case and the power-law case will produce, respectively, either flat or decreasing energy 
spectra [0, 

Having said this we will focus our attention to the case where a pure de Sitter phase is 
followed by a stiff phase and we will comment, where appropriate, on the other possible cases. 
Therefore the model we want to investigate is given (in conformal time) essentially by three 
phases: 



(l + /5)r/i + r7]/3 



The subscript i, s and r simply stand for inflationary, stiff and radiation dominated phases. An 
important feature of Eq. (|2.6| ) is the continuity of the scale factors (and of their first conformal 
time derivatives) in the matching points rji and 7]^. Notice also that the generic exponent f3 
specifying the dynamics during the stiff phase is trivially related to the, previously introduced 
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/3=^ = ^. (2.7) 
^ 1-a 37 + 1 ^ ^ 



a and 7 parameters, namely 



Of course the three phases specified by Eq. ( |2.6| ) are usually complemented by the transition 
to matter- domination occurring at rjdec- For 7] > rjdec the scale factor evolves parabolically (i.e. 
O'm.iv) ~ v'^) ill the matter epoch. Our main goal in this paper will be to study the gravitational 
wave spectra in these models and to study which kind of constraints arise by changing Hi and 
7. Notice that the time of the infiation-stiff phase transition is rji = [ai{ti)Hi{ti)]^^ . 



3 Gravitational wave spectra from stiff phases 

The evolution of the scalar, vector and tensor fiuctuations of a given background geometry 
can be directly discussed by perturbing (to second order in the amplitude of the fiuctuations) 
the Einstein-Hilbert action. An important property of the metric perturbations in FRW back- 
grounds of the type defined in Eq. (|2.1D is that scalar, vector and tensor modes are decoupled 



p4| (this feature does not hold in the case of anisotropic background geometries ||23|). This 
fact means that, by defining the fiuctuations of the background metric as 

9f.u{x, rf) = g^^{ri) + 5^^^(x, r/), (3.1) 

the metric fiuctuation can be formally written as 

5g^,{x, T]) = 5g%\x, r]) + 5g%\x, r]) + 5g%\x, 77), (3.2) 

where S", V, T stand respectively for, scalar, vector and tensor modes. This classification refers 
to the way in which the fields from which 5g^j^p{x,ri) are constructed, change under three di- 
mensional (spatial) coordinate transformations on the constant (conformal)-time hypersurface. 
Since 5g^i,{x,ri) (being a symmetric four- dimensional tensor of rank 2) has ten independent 
component the scalar, vector and tensor modes will be parametrized by ten independent space- 
time functions. More specifically, scalar perturbations will be parametrized by four indepen- 
dent scalar functions, vector perturbations by two divergence-less three-dimensional vectors 
(equivalent to four independent functions). Pure tensor modes of the metric (corresponding to 
physical gravitational waves propagating in a homogeneous and isotropic background) can be 
constructed using a symmetric three tensor hij satisfying the constraints 

V = 0, /ii = 0, V^7i,, = (3.3) 
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( Vj denotes the covariant derivative with respect to the three-dimensional metric). Notice 
that, with our conventions, Sg'jfJ = h^y and 5g^^^ = —h^". The hne element perturbed by the 
tensor modes can then be written as 



(3.4) 



(where 'jij is the spatial background metric). 

Summarizing, we have four functions for scalars, four functions for vectors and two func- 
tions (the two polarizations of hij) for the tensors. In total there are ten independent degrees 
of freedom describing the fluctuations as required by the tensor properties of the original (un- 
perturbed) metric. Moreover the conditions expressed in Eq. p.3| ) imply that the two physical 
(independent) polarizations of do not contain any pieces which transform as scalars or 
vectors under three-dimensional rotations. As a consequence of its definition, hij is directly 
invariant under infinitesimal coordinate transformations preserving the tensorial character of 
the fluctuations [§ 

In order to obtain the evolution equation of the metric fluctuations there are at least two 
different procedures. First of all one could think to perturb (to first order in the metric fluctu- 
ations) the Einstein equations. On the other hand one could also perturb the Einstein-Hilbert 
action 



S = - Jtt I d'^X\f-9R, R = 5'"^^a/3, 9 = det 



6Pp 



Ip = Mp\ (3.5) 



to second order in the amplitude of the metric fluctuations ||25|| . 

It is convenient to notice that the two procedure are certainly equivalent but the perturba- 
tion of the action provides more informations since it allows to isolate (up to total derivative 
terms) the normal modes of oscillation of the system which one might want to normalize to the 
value of the quantum mechanical fluctuations. By perturbing the action given in Eq. ( p. 51) to 
second order in the amplitude of the tensor fluctuations we obtain that, up to total derivatives 



'''''''' - 2k I 



dahijdph'^f 



v/3 



(3.6) 



(we remind that in this and in the following formulas the shift from upper to lower spatial 
indices [and viceversa] is done by using the spatial background metric and its inverse). 
For a wave moving in the = z direction in our coordinates one has = h\ = —h\ and 
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= h\ = h\ and the perturbed action becomes: 



5(2) 5{T) 



12/2, 



j c/^xa/^ dahQ)df3h(sg"'^ + dahf^dph^g' 
By now varying the perturbed action we get the evolution equation for each polarization: 

h@ + - V^/ie = 0, + 3Hh^ - W^h^ = 0. (3. 



(3.7) 



From Eq. (|3.7|) it is also possible to deduce the form of the canonical normal modes, namely 
those modes whose Lagrangian reduces to the Lagrangian of two minimally coupled scalar fields 
(in flat space) with time dependent mass terms. Defining 



we can write, from Eq. (|3.7| ), the action for the canonical normal modes 
where (always up to total derivatives) 



1 



(3.9) 



(3.10) 



(3.11) 



iVai3 = diag(l, — 1, — 1, — 1) is the fiat space-time metric and H = [log a]'). From Eq. (|3.11|) 
we can derive the evolution equations of the canonical normal modes by taking the functional 
variation with respect to fi,^ and /x© with the result that: 



0. 



(3.12) 



f^l ~ vV® — = 0, _ 

a a 
The normal modes obtained in Eq. (|3.19| ) can be now canonically quantized. First of all we 
define the canonical momenta 

dL dL , . 

leading to the Hamiltonian 



Q 



(3.14) 



We then impose the (equal time) commutation relations between the corresponding field oper- 
ators: 



i5'\x-y), 



iS^{x — y). 



(3.15) 
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We want to remind that the way we got the Hamiltonian is quite naive. In fact there is a 
(formally more correct) way of getting the Hamiltonian for the perturbations ||2^. One should 
in fact start with the Hamiltonian in the superspace, fix the background, and then obtain the 
Hamiltonian for the perturbations together with four (independent) constraints (two for scalar 
perturbations and two for vector perturbations). Therefore, in a fully consistent Hamiltonian 
approach to perturbations there are no constraints arising from the tensor modes. This is the 
reason why our naive approach leads to the same result of the Hamiltonian formalism for what 
concerns the tensor modes. In the case od scalar fluctuations, however, our approach would 
be less correct since, in the scalar case, the off-diagonal components of the Einstein equations 
receive contribution leading to constraints. In order to get the Hamiltonian for the scalar 
modes (and the related algebra of the constraints) the (Hamiltonian) superspace approach is 
more compelling. We point out, in any case, that also in our present approach it is possible to 
get the correct result for the Hamiltonian of the scalar modes by inserting, in the perturbed 
action of the scalar fluctuations, the constraint equation arising from the {Oi) components of 



the Einstein equations |27 



Promoting the classical Hamiltonian quantum mechanical operator the evolution equations 
for the field operators become, in the Heisenberg representation. 



i/i^ — 


/■ ■ Q 










1 '^'^(g) ~ 


Q 



(3.16) 

Notice that (as it has to be) the evolution equations in the Heisenberg representation are exactly 
identical (for the field operators) to the ones derived for the classical fields in Eq. (|3.12| ) once the 
explicit expressions of the conjugated momenta (i.e. vr^ = fi'^ and tt^ = fi'^) are inserted back 
into Eq. ( |3.16| ). We also point out that the hermitian Hamiltonian of Eq. ( |3.14| ) is quadratic 
in the field operators and it belongs to a general class of time- dependent Hamiltonians widely 
used in quantum optics ^ in the context of the parametric amplification of the vacuum 
fluctuations of the electromagnetic field through laser beams. The same discussion we performed 
in the Heisenberg picture can be easily translated to the Schroedinger picture where the final 
state of the evolution is a many particle state unitarily connected to the initial vacuum (a 
so called squeezed state |^). Having fixed these standard notations we can expand the 
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operators in Fourier integrals 

1 



(27r)3/2 

1 

(27r)3/2 



-ik-x 



-ik-x 



(3.17) 



(recall that, in our notations, h^^^{x,T] = ^/6[lp / a]jl(£i^{x, r])) . From ( p.l5| ) a© and obey the 
following commutation relations 



a®(A;),a^(A^') 



5^^\k-k') 



d(x,{k),al,{k') 



5^^\k-k') 



d(s{k),d(^{k'] 



0. (3.18) 



Therefore the evolution equation for the Fourier amplitudes ^^{k.rf) and fi(^{k,ri) become 



a 



(3.19) 



(where, for the Fourier amplitudes, we define fi{k,rj) = ah{k,rj)). 

By taking the functional derivative of the action reported in Eq. (|3.6|) with respect to the 
background metric we indeed find an effective energy-momentum tensor of the fluctuations 
which reads: 



T, 



(3.20) 



If we define the vacuum state vector |0©0g,) which is annihilated by a® and d^ we obtain that 
the energy density of the produced gravitons will be given by 



1 



d-'k 



h'Jk,r^)\' + Kik,r^)\' + k'lh^ik^v)]' + Kik,v)\ 



(3.21) 



In order to compute the GW spectra produced by the transition of the background from an 
inflationary phase to a decelerated, stiff, phase we have to solve the evolution equation ( p.l9| ). 
From Eq. ( 3.19 ) we clearly see that the evolution equation of the (tensor) normal modes of the 
geometry is determined not only by but also by H. Therefore, depending upon the sign 
of H different inflationary models will give different (large scales) spectral distributions of the 
amplified gravitational waves. 
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f(ri)=a2(H + 2H ) 



Figure 1: The "effective potential" a" /a appearing in the Schroedinger-hke equation ( p.l9| ) is 
reported as a function of the conformal time coordinate. In the standard scenario, the inter- 
mediate (stiff) phase is absent and, therefore, all the modes k will reenter during the radiation 
epoch leading (at high frequencies) to the usual Harrison- Zeldovich (flat) energy spectrum. In 
our case the modes re-entering during the radiation dominated era {k < rj^^) will always have 
flat energy spectrum. At the same time, some modes going under the potential barrier during 
the de Sitter phase (z-stage) and re-entering the horizon during the stiff phase (s-stage) will lead 
to a modification of the high frequency branch of the spectrum which will become mildly increas- 
ing (see Eqs. (|3:30|)-(|33l|) and Eq. (|33^)). We stress that at large (lO-^^Hz <u < IQ-^^Hz) 
and intermediate (10~^^Hz < u; < uj,,.) (present) physical frequencies, the spectrum does not 
change with respect to the usual results [0, |^, ||, |^ since the corresponding modes will not 
"feel" the presence of the stiff phase. 



In the three phases defined in Eq. (plGp the evolution equation ( p.l9|) reads : 



fi = 0, 



/i = 0. 



7]i <r] < r]r, 



fx" + k'^jj, = 0, 



Vdec <r] <r]r 



(3.22) 



(notice that we dropped the subscript referring to each polarization). The solution of the 



evolution equations in the three regions can be written in terms of Hankel functions |34 



Since we are dealing with the Fourier amplitudes of the normal modes of oscillation we normalize 
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them directly to the quantum mechanical noise level (~ l/\/k) for rj < —rji: 

1 

b+Hl^\v) + b^Hl'\v 



'2k 



Vdec <r] <rir, (3.23) 



where 



2z/=|2/3-l|, u = kr], v = k[r] + {p + l)r]i]. (3.24) 



Notice that the Wronskian normalization of the solutions imposes that — = 1 and 
|c+p — |c_p = 1. See also Fig. |l| for a pictorial description of the amplification process. By 
computing the coefficients 6_ and c_ we can have an estimate of the energy spectrum of the 
produced gravitons re-entering the horizon, respectively, after rji (i.e. right after the completion 
of the inflationary phase) and after 1]^ (i.e. immediately after the beginning of the radiation 
dominated epoch). Notice that rir appearing in Eq. (|3.23|) is not free parameter since the 
duration of the stiff phase ending in rjr is determined essentially by the backreaction of the 
modes re-entered right after t]i. In order to further clarify this point we have to compute the 
GW energy spectra produced right after rji. The energy density of the produced gravitons per 
logarithmic interval of longitudinal momentum can be obtained from Eq. ( p. 211 ) : 

dioguj 

where uj = k/a. Sometimes p{uj,t) is called logarithmic energy spectrum. Notice that, one can 
deflne the energy density of the gravitational wave background (in a semiclassical description) 
through the energy momentum pseudo-tensor of the amplifled tensor modes. In the approach of 
the energy momentum pseudo-tensor the quantum mechanical expectation values are replaced 
by ensemble averages over a distribution of stochastic variables whose two-point function has 
to be specified separately by requiring that each Fourier amplitude of the tensor modes is 
independent on the others in momentum space. This "stochastic" condition on the Fourier 
amplitudes is the result, in a fully quantum mechanical approach, of the unitarity of the graviton 
production process driven by the hermitian Hamiltonian (|3.14|) . We stress that, if our initial 
state is not the vacuum (but for example a non-pure state characterised by a specific density 
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matrix 
different. 



the expression given in Eq. (|3.21|) for the averaged energy density will also be 



From Eq. ( p.21| ) and using the definition (|3.25|) we have that in each phase of the model 



1 



the energy spectra can be written as: 

p{uj,t) 



27r2' 
1 





,4 

) \l>-? 


"Oi" 








. a _ 






'ai' 












. a _ 



V > Tjr. 



(3.26) 



In order to compute the coefficients and |c_| we use the sudden approximation, namely 
we will mainly consider the amplification of the modes kr]i<l A possible technique is to 
match, in rji and rjr the exact solutions given by Eq. ( |3.23|) and their first derivatives. Therefore, 
imposing the conditions 



PniVr) = PlIl{Ur), p'lliVr) = p'jjjiUr), 



(3.27) 



we obtain an expression for the amplification coefficients which we evaluate in the sudden 
approximation: 
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UJ 
UJr 



-I- 



'2 f UJ 



UJi 



Ur < UJ < Ui, 



UJdec < UJ < UJr 



(3.28) 



(notice that Ui^r = ^'7i,r and Vi^r = ^[^i,r + (/? + l)^i])- These expression are correct as long as 
the effective equation of state parameterizing the evolution of the sources in the intermediate 
phase has 7 < 1. For 7 = 1 (corresponding to jS = 



a 



I and = 0), the solution of 
Eq. ( 3.22| ) in the intermediate phase (i.e. what we call pii{v) in Eq. ( |3.23| )) is given in terms 



of Hq ' (f), and, therefore, the amplification coefficients get logarithmically corrected (in the 
sudden approximation), and the final result is : 

l^-l - ( — ) log 



UJ 
UJl 



UJi 

UJr 



UJr < UJ < UJl, 



UJr 
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c<;\ 2 f UJ \ 2 ( UJl' 
— — log — ;> 



UJl J 



UJr 



UJdec < UJ < UJr 



(3.29) 
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Notice that uJdec — 10~^^ Hz is the frequency corresponding to rjdec-, namely to the moment of 
the transition to the matter-dominated epoch. Of course the gravitational wave spectra will 
have a further (fourth) branch corresponding to the modes re-entering the horizon during the 
radiation dominated epoch This infrared branch of the spectrum concerns frequencies 

ujq < UJ < UJdec where ujq ~ 10"^*^ Hz is the frequency corresponding to the present horizon. It is 
well known that in the infra-red branch the spectrum turns out to be steeply decreasing as a;~^ 

0. Since the purpose of our investigation are the effects of the produced gravitons during 
the stiff phase we will neglect, in the following considerations, the infrared branch. 

By inserting the obtained expressions of the amplification coefficients in the expression for 
the energy spectrum of Eq. ( p.26|) we get the main result of this Section 

1 
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Again, this expression is logarithmically corrected in the 7 
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UJr < UJ < UJi, 
UJdec < UJ < UJr- 

1 limit: 

UJr < UJ < UJi., 
UJdec < UJ < UJr. 



(3.30) 



(3.31) 



Concerning Eqs. ( |3.3(J| ) and ( |3.31| ) few comments are in order. First of all we can notice that 
the energy-spectra are scale-invariant for uj < uJr This branch of the spectrum corresponds 
to modes which went out of the horizon during the initial de Sitter phase and re-entered 
during the radiation dominated phase. There is, therefore, no surprise for this behaviour which 
does correspond to the usual Harrison-Zeldovich spectrum of the stochastic GW backgrounds. 
Moreover, we point out that for uJr — > uji (i.e. point-like stiff phase) we reproduce the usual 
results well known in the context of the theory of graviton production when a de Sitter phase 
is suddenly followed by a radiation dominated phase [|12], ^ For the modes leaving the 
horizon during the de Sitter phase and re-entering during the stiff phase {uJr < uj < uoi) we 
get the curious result that the energy spectrum increases as a function of the frequency. This 
feature of our result can be easily understood by bearing in mind that the ultraviolet branch 
of the spectrum goes as 

3(1-7) 
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(3.32) 
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If I < 7 ^ 1 we have that 2z/ < 1 and the spectrum is always mildly increasing with maximal 
slope {uo/uji) (for 7 = 1). The minimal slope corresponds to the case 7 = 1/3 (flat case) where 
the fiat spectrum is recovered for all the (present) frequency range. This peculiar behaviour 
is indeed not so strange. This peculiar feature of stiff models was indeed noticed long ago 
2U|. In our case the only difference is that prior to the stiff phase there is a de Sitter phase 
and, therefore, the calculation of the amplification coefficients involves a further transition. 
This further transition modifies the high energy behaviour of the spectrum which still increases 
but more mildly if compared to the case where the de Sitter phase was absent |2^. If the 
inflationary phase is not de Sitter but power-law (i.e. < 0) or superinfiationary [H > 0) 
the spectra of the modes re-entering in the radiation epoch will be crucially modified and their 
amplitude will be always subjected to the large scale constraint {Hi<10~^Mp) provided the 
inflationary phase is either de Sitter-like or power law. In the case of superinflation (leading 
to increasing energy spectra), the nucleosynthesis constraint will always be the most stringent 
Sy, ^ since it involves the integrated energy density and it applies at all the frequencies. 



one 



The background energy decreases, during the stiff phase, as HfMp{ai/a)^^'^^^\ whereas 
the radiation stored in GW decreases like (ai/a)^ and therefore, at some stage, the graviton 
radiation will become dominant. To compute precisely this moment we have to integrate the 
graviton spectrum over all the modes, insert it back in the Einstein equations and, finally solve 
the modified Einstein equations. By defining e = Hi/Mp, a simple argument based on the 
calculations reported in the previous Section shows that, since the energy spectra are increasing, 
the most significant contribution of the hard gravitons to the energy density (integrated over 
the whole spectrum) will occur for ~ cji. This energy density will turn the stiff background 
into radiation at a critical value of the scale factor and the stiff fluid will 

correspondingly turn radiation (i.e. 7 1/3). Taking now into account that during the stiff 
phase a{t) ~ t'^/i^i'y+'^)] have that, according to our estimate the back reaction effects will 
become signiflcant at a curvature scale Hr ~ e^^^'^'^^^^^"'~^^ Mp. For example, if we take 7 = 1, 
the curvature scale at which the transition to radiation takes place is H,. ~ e^Mp. In frequency, 
the length of the stiff phase (with 7 = 1) will then be Ur/uJi ~ rjiai/rj^ar ~ e^. Thus, if one 
wants to reach (within a stiff phase with 7=1) the hadronic curvature scale Hhad ~ 10~^° Mp 
we should flne-tune e ~ 10~^°. 

The maximal scale Hi is certainly constrained. If a de Sitter phase is immediately followed 
by a radiation dominated phase, the GW energy spectrum is flat for uj^ec < oj < uJi (notice 
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that, today, Udec ~ lO^"*^^ Hz and ui ~ lO^-'^v^ Hz). For uq < u < ojdec the spectrum decreases 
and, therefore, e<10^^ which is exactly the bound reported in Eq. ( p. . 1|) coming from the tensor 
contribution to the CMBR anisotropy. 

If the de Sitter phase is followed by a stiff phase the spectra are growing for ojr < oj < uJi. 
Take, for instance, the case 7 = 1 where p{uj,t) ~ Hf{uj/uji) up to logarithmic corrections. 
In this branch, since the spectrum increases, the most significant constraint comes from the 



bound energy density in relativistic degrees of freedom at the nucleosynthesis epoch pOl 31 



Thus the bounds on e might be a bit different (even if not by much due to the very mild 
increase of the spectral energy density). Today, in the range Udec < uj < ujr the spectral energy 
density (in critical units) would be flcwi^jt) ~ 10~^e^{ujr/u!i). Now the large scale bound 
imposes that, in this phase, flcw^^^^^^ ■ Taking now into account that Ur/uJi ~ (for the case 
7 = 1) we have that flcwi^yt) ~ e^- Imposing now e< 10^^ we get that, for Udec < uj < Ur, 
^Gw{uj,t) < 10^^^. The argument we just discussed can be easily extended to the entire class 
of stiff models (i.e. 1/3 < 7 < 1). 

In conclusion the possibility of having a stiff, post- inflationary phase implies necessarily a 
"reheating" driven by gravitational waves re-entering the horizon during the stiff phase. In 
order to have a significantly long stiff phase, however, fine-tuning is strictly required making 
these models perhaps less attractive. We will elaborate on this point in the next Section. 



4 Back reaction effects 



In order to estimate the length of the transition between the stiff phase and the radiation 
dominated phase induced by the hard gravitons, we re-write Eq. ( |3.19| ) in a slightly different 
form, namely 



dt 



where fife(r) 
as 

with 



/iffi «?, = 0, dri = — = a dr, (4.1) 
-guj = ka? and h = a~^fi. A formal solution to this equation can be written 

1 



h{k,T) 



h±{k, t) = exp 



C+h4k,T) +C^h4k,T) 



±i / Qudr 



(4.2) 
(4.3) 
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(notice that C+ and C_ are complex functions of r; thanks to the Wronskian normahzation 
condition we also have that |C+p— = !)• By inserting the solution given in Eq. ( [4.2|) back 
into Eq. ( [4.1|) we can obtain an evolution equation for the two, time-dependent coefficients, C± 
161 

dC+ 1 dlogfii. , r, _ dC^ 



2 dr + " 



l rflogl]fc ^2 

2 dr ' ^ 



(4.4) 



dr 2 dr dr 
(from now on we will drop the subscript referring to the two polarizations). We are now going 
to solve these equations for the modes krii<l. It is easy to show that this sudden approximation 
used in the previous section corresponds, in the language of Eq. ( [4. 1] ) to the small J'^^ (irfifc(r) 
limit. In fact 

/ Vtk{r)dr = / ka dr = / kdr] ~ kr]i (4.5) 

(the last two equalities follow from the definition of dr in terms of the conformal time coordinate 
: dr = drj/a^). 

In this approximation we can expand the h± appearing in Eq. 
order in / drQ^, 



and we find, to first 



dC+ 
dr 



1 d log Qk 

2 d^ 



l + 2i Qkdr 



dC^ 
dr 



1 d log Vlk 
2 



dr 



l-2i Q.dr 



(4.6) 



By now linearly combining two previous equations we find: 



d 

dr 
d 

dr 



1 d log Qk 

2 dr 
Idlogfi 

~2 rfT" 



+ 



+ ... 



(4.7) 



In Eq. ( [4.7[ ) the ellipses stand for other terms which are of higher order in /^^ VLk{r)dr ~ krji 
and which are negligible in the sudden approximation. The solution to Eq. ( [4. 71) can be easily 
found in terms of two (arbitrary) complex coefficients 

1 



C+ + C_ = 2QiA/l^fc, C+ - C_ = 2Q2 



/a 



As we previously mentioned, the Wronskian condition imposes that IC4 



(4i 



1. Now 



this last condition has to hold order by order in krji and, therefore, inserting the C± of Eq. 



into \C^ 

C4 



|C_p = 1 we obtain a condition on Qi and Q2 valid to first oder in krji. 



Q2 



Q2 



QIQ2 + QiQl 



(4.9) 
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Since the amount of amplification of the GW is essentially determined, in this approach, by 
we fix Qi and Q2 by requiring that at the time 771 (or ti) \C^{k,ti)\ = 0. Thus, from Eq. 
we get that 



IQ2 



Using Eq. 

PGwit) 



2V^fc(ti) 

^ and summing the polarizations we get that 



(4.10) 



d ujijj 



|C+(cu,t)p + |C_(cu,t)p 



1 



ijjd uj 



2|C_(cj,t)|2 + l 



(4.11) 



(where in the last equality we used the Wronskian normalization condition). To the lowest 
order in kr]\, using Eqs. ( [4. 91 ) and (|4.8|), and expressing all the quantities in terms of the 
corresponding physical momenta a;, we get 



1 



16 7r3 



+ 



(4.12) 



Now we want to compute how much energy is present at a generic time t inside the horizon 
during the stiff phase. First of all we will outline the formal solution of the problem and 
secondly we will do an explicit calculation. 

The total energy density in GW at a generic time t will be given by the explicit momentum 
integral indicated in Eq. ( |4.12| ). This is however not the end of the story. At any given time 
different GW will re-enter during the stiff phase and therefore the total energy stored in GW 
which should appear at the right hand side of the Einstein equations is the sum over all waves 
re-entering at different times after ti. This means that the quantity we we should compute is 
[p2| the total energy density of the gravitational waves re-entered at some generic time t during 
the stiff phase. We then have that the total energy density of the graviton background is given 
by the modes re-entering at t summed to the (red-shifted) energy density of those modes which 
re-entered during the period ti < f < t: 



Ptotit) 



a(t') 



a{t) 



duj„ 

a~dt' 



-dt' 



(4.13) 



Concerning Eq. ( [4.13|) few comments are in order. First of all PGwiy^mit) denotes the energy 
density in GW (see Eq. ( |4.12|) ) integrated until a "running" ultra-violet cut-off uj^nif) ~ j ~ 
H{t). The partial derivatives appearing in Eq. ( ^.13| ) simply reflect the fact that we are 
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summing up the energy of GW re-entering the horizon at different moments t > ti. In order to 
do this we have exactly to shce the horizon in many infinitesimal portions, and compute, for 
each portion, the corresponding increment in the GW energy. In computing the increment in 
the GW energy with respect to the cut-off we keep the scale factor constant (and this is the 
reason of the subscript appearing in Eq. ( 4.13| ). We want now to compute the frequency integral 
of Eq. (|4.12 ). As it is well known this integral is divergent in the limit of large frequencies. 
This is simply a consequence of the fact that the vacuum leads to an infinite energy density 
which has to be properly subtracted. The spectrum of the vacuum fluctuations can be simply 
obtained by putting \C-{uj,ti)\ = in Eq. ( [4.11| ). This limit corresponds to the absence of 
amplification, and therefore the only fluctuations contributing to the energy density are the 
ones associated with the vacuum modes with logarithmic energy spectrum proportional to u'^. 
In principle, in our case we have a (physical) ultraviolet cut-off in the spectrum provided by the 
scale where inflation stops. This maximal frequency is ^ We then expect our results to 
be insensitive to the particular renormalization scheme. Since however we want to have (in Eq. 
(|4.13|) ) the possibility of a cut-off running with time we will examine briefly this issue which 
was actually investigated in the past, for homogeneous cosmological backgrounds, within (at 



least) two different approaches. In Ref. [^] this problem was tackled using a regularization 
scheme strongly reminiscent of the Pauli-Villars method. In Ref. [^] the same problem was 
discussed within the so called adiabatic regularization scheme (see also P2[). The two methods 



were shown to produce equivalent results fS^]- The purpose of this paper is not to check which 
regularization scheme is better to use in curved space. Our approach is more pragmatic: we 
want to get an estimate not only of the maximal duration of the stiff phase but also of the 
transition time between the stiff and the radiation dominated epoch. We want to know how 
long will it take for the stiff fluid to turn into radiation. In this spirit we will firstly of all use the 
adiabatic regularization scheme (without including higher order subtractions involving double 
time derivatives of the Hubble parameter). Secondly we will re-compute our back reaction 
effects on the background evolution without including the subtractions and using a (naive) 
cut-off regularization. 

In the adiabatic scheme the regularized energy density reads 

1 



Pewit) 



(4.14) 
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where 



+ 



2\uj 



(4.15) 



We notice that Fi{ijj,t) is simply the non-regularized contribution and it is exactly equal to 
the one we computed Eq. ( [4.11| ). On the other hand, F2{uj,t) comes from the adiabatic 
regularization and it does contain the divergent contribution which we want to subtract. 
Notice that, in the expansion there are terms with higher derivatives since 



+ 2m - 2HH - AH^H 



(4.16) 



Now, our underlying theory is the Einstein theory (with only linear curvature terms in the action 
( |3.5| )). Consequently, it is not compatible to include subtractions involving four derivatives. 
Therefore we will not include them in the subtraction. We can integrate the energy density 
keeping a running ultra-violet cut-off ~ H{t)<Ht^ 



PGwit) 
with the result that 
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(4.17) 



(4.18) 



We are now ready to include the effect of the produced gravitational waves in the Einstein 
equations. The evolution of the Hubble parameter will be given (in the conformally flat case) 
by the following integro-differential equation: 



MlH^ - p,{t) = - f 
Jti 

Ps{t) = HlMl 



dp, 



GW 




GW 



doJn 



-dt' 



(4.19) 



(4.20) 



In Eq. ( |4.20| ) ps is just the energy density of the stiff background which decreases faster that 



for any 7 > 1/3. Notice that the right hand side of Eq. (|4.19|) is noting but the first of 



the FRW reported in Eq. ( p. 21) . In the absence of graviton creation the left hand side of Eq. 
(|4.19|) would be just zero, whereas in the presence of graviton re-entering the horizon at any 
time after ti this second term receives a non vanishing contribution. 
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An useful way of re-writing the integral appearing in Eq. ( |4.19| ) is by changing integration 
variable from t' to H' = H{t'), and in this way the total energy of produced gravitons becomes: 



Ptot(i) 



4vr2 
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l4 
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ai 
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H' 



>dH' 



(4.21) 



Notice that dum/dH ~ 1 since ujm ~ H when the given mode crosses the horizon. 

It is clear that Eq. ( [4. 201 ) represents a complicated integro-differential equation which 
cannot be solved exactly. One possible way of dealing with this problem is to transform it into 
an ordinary differential equation by solving the integral using the scale factor of the stiff phase. 
In other words the integral at the right hand side of Eq. ( [4.19|) can be viewed as a perturbation 
to the solution of the FRW equations ( |2.2| ) in the absence of graviton creation. Using this 
procedure we insert the stiff scale factor into Eq. ( [4.20| ) and we compute the effect of the 
gravitons re-entering after ti on the background metric. Thus using [a{t') / a{t'^)] ~ ^' /ti\^/^^'^'^^\ 
ijjm ^ ^/t' and H' = H(t') ~ 2/[3(7 + l)t'] into Eq. ( [4.21|) , our integro-differential equation 
becomes 



dyj 



9(7 + 1)2 

2(1+37)^ 
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/i(7) 



A(^)= /i(7) 



1 2g + 9(l + 7) 
247r2 (7 + l)(37 + l)2' 
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(4.22) 



(4.23) 



(recall that e = Hi/Mp ). In Eq. ( [4.23| ) there is also the parameter q which needs to be 
explained. As we said before we regularized the energy density using the adiabatic regularization 
scheme. Now, if we set g = 1 we automatically include in Eq. (|4.23| ) the subtractions coming 
from the adiabatic regularization. If we set g = we practically use for the calculation of the 
energy density of the gravitational waves the non- regularized energy density. Thus', if g = 1, 
Pgw is used in Eq. ( [4.201 ). If g = 0, Eq. ( |4.2CI| ) is computed by using pcw which does not 
include the subtractions of the divergent terms. As we discussed previously to set g = 1 or 
g = does not change the numerical solution which we are going to describe. It should be 
actually be borne in mind that we have a physical cut-off provided by the class of models we 



are discussing and which is set by ki 



The effect of the subtractions (encoded in the 
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Figure 2: We plot z'{y) in the case e = 10^^. Recall that z = log ^ and y = log ^. Therefore, 
we have that a{y) = z'{y). We plot the solutions for four different values of the initial 7 
parameter. We integrated Eq. ([4.23|) from t = ti (i.e.?/ = 0) up to y = 120. Take, for instance, 
the case 7 = 1 corresponding to a{ti) = 1/3. We see that, thanks to the gravitational waves 
contribution the background undergoes a phase transition towards a radiation dominated epoch 
with z'{y) ~ o{y) ~ 1/2. In the full lines is reported the case g = (no subtractions of the 
ultra-violet divergences) whereas in the full lines is illustrated the case q = 1 includeding the 
subtractions. 
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Figure 3: The time evolution of the equation of state is reported. Recall, for comparison 
between the two pictures that 7(?/) = 2/[3z'(y)] — 1 and that a{y) = z'{y). The four curves 
do correspond to the four cases already discussed in Fig. We notice that starting with stiff 
equations of state at t = ti (i.e., from top to bottom, 7(0) = 1, 5/6, 3/4, 2/3) we are attracted 
towards 7 = 1/3 for large y. 



choice of q) is illustrated in Fig. ^ where the solution of Eq. ( f4.23| ) is reported as a function 
oi y = logt/ti. We can see that to include the subtractions affects the transition regime (when 
the gravitational radiation starts to dominate the background) but not the asymptotic regime. 
Notice that from the effective evolution of a{y) we can derive also the effective evolution of 
7(2/). In Fig. ^ the cosmic time evolution of 7(2/) is illustrated for the four different cases of 
7(0) discussed in Fig. ^. The back reaction which forces the equation of state to pass from 
its original stiff value (for y — > 0) to 1/3 typical of a radiation dominated (perfect) relativistic 
fluid. 

It is also clear that, by lowering e the transition to the radiation dominated epoch might 
be delayed (paying, of course, the price of a fine-tuning). In Fig ^ we report the solution of 
Eq. ( [4.23|) for different values of e in the case 7 = 1. As usual we can either see the transition 
either in terms of 7 or in terms of a. 

Looking at Fig. ^ and ^ we can also see that the transition to the radiation dominated 
phase does not occur instantaneously but it takes place in a finite amount of time which turns 
out to be quite substantial (of the order of 20 time e- folding). 
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Figure 4: The solution of Eq. ( ^.231 ) is reported for the case 7 = 1. If e = 10"'* (as required by 
BBN constraints) we have that the equation of state of the fluid driving the expansion starts to 
deviate significantly from the stiff one already for ?/ ~ 4 x 10^ corresponding to Hr ~ lO^^^Mp 
(notice that always used neperian logarithms). 



5 Discussion and Conclusions 

In this paper we investigated the possibility of stiff epochs occurring immediately after an 
inflationary phase. We found that, if these phases ever existed, they are constrained by the 
production of GW. We computed the associated energy spectra of the produced gravitons 
and we also discussed the associated back reaction effects. For what concerns specifically the 
theoretical implications of GW backgrounds of stiff origin, we can say that the possibility of 
having blue spectra at high frequencies looks certainly promising. In this class of models "blue" 
spectra arise quite naturally. If the maximal inflationary curvature scale is taken to be of the 
order of 10~^ Mp we can also see from our results that the amplitude of the GW (logarithmic) 
energy spectrum will be larger, at high frequencies, than the inflationary prediction obtained 
in the absence of stiff phases. The backreaction effect will make the Universe dominated by 
radiation and the corresponding "graviton reheating" energy scale in the range 1-10^ TeV. 

The duration of the stiff phase can be long only if the maximal scale where inflation occurs 
is fine-tuned to be much smaller than 10~^Mp. We also found that the transition regime (where 
the stiff equation of state with 7 > 1/3) is replaced by 7 = 1/3 is quite long. Our analysis has 
certainly different limitations. We mainly focused our attention on the back reaction effects 



associated with high frequency gravitons which behave effectively like radiation [2C, 21, 22 . 
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The second limitation of our analysis is that we did not discuss the amplification of the scalar 
fiuctuations. We did not do this for the simple reason that the scalar fiuctuations are much 
more sensitive to the particular dynamical model used in order to implement a stiff phase. 

A conservative conclusion of our investigation is that the the duration of a stiff, post- 
inflationary phase is crucially determined by the energy density of the inhomogeneity which 
were amplified during the inflationary phase and re-entered in the stiff phase. In the framework 
of a particular model the back reaction effects of tensor (and scalar) fluctuations should be 
analyzed. On one hand the produced inhomogeneities could offer an original mechanism for 
re-heating the Universe. On the other hand they can forbid an arbitrary duration of the stiff 
phase. 
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